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A c-axis magnetotransport and resistance noise study in Lai.97Sro.o3Cu04 reveals clear signatures 
of glassiness, such as hysteresis, memory, and slow, correlated dynamics, but only at temperatures 
(T) well below the spin glass transition temperature T sg . The results strongly suggest the emergence 
of charge glassiness, or dynamic charge ordering, as a result of Coulomb interactions. 
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The role of heterogeneities observed in most hole- 
doped high-temperature superconductors (HTS) is one of 
the major open issues in the field [H, 0] . In weakly doped 
Mott insulators, such as HTS, charge heterogeneities are 
expected to arise due to the existence of several compet- 
ing ground states Q, and may be even to exhibit glassy 
dynamics 0). Even a small amount of disorder may fa- 
vor glassiness over various static charge-ordered states 
@| . Experiments in hole-doped HTS suggest @, 0, S 
that glassiness of both spins and charges emerges with 
the first added holes and evolves with doping x. While 
spin glass (SG) behavior is well established at low T, 
the evidence for glassy freezing of charges is not conclu- 
sive. Hence, alternative, bulk probes of charge dynamics 
are needed to explore the nature of the ground state. 
We present a novel study of the charge dynamics in a 
lightly doped I^-xSr^CuO^ (LSCO) using a combina- 
tion of transport and noise spectroscopy that proved to 
be a powerful probe of dynamics in other glasses ToL 
We find several clear signatures of glassiness at T <C T sg . 
The data strongly suggest that the doped holes form a 
dynamically ordered, cluster glass as a result of Coulomb 
interactions. 

In LSCO, the prototypical cuprate HTS, the three- 
dimensional (3D) long range antiferromagnetic (AF) or- 
der of the parent compound is destroyed above x « 0.02, 
but 2D short range AF correlations persist [T3| • In partic- 
ular, as a result of hole doping, Cu02 (ab) planes develop 
a pattern of AF domains that are separated by antiphase 
boundaries 13|, [bj, [l5| ■ Since the Dzyaloshinskii-Moriya 



interaction induces slight canting of the spins in Cu02 
planes towards the c axis, there is a weak ferromagnetic 
(FM) moment in the be plane associated with each AF 
domain, such that the direction of the FM moment is 
uniquely linked to the phase of the AF order [3, flij ]. 
The interplane exchange favors staggered ordering of 
those FM moments in the c direction. At low enough 
T < T sg (x), the system freezes into a SG [H, E3, El that 
extends into the superconducting (SC) phase for x > 0.05 



[171 ] up to optimal doping [18| [Fig.QJa)]. Various exper- 
iments on lightly doped LSCO (e.g. Refs. [1,0]), includ- 
ing transport studies (lij , were interpreted in terms of the 
hole-poor AF domains separated by the hole-rich regions 
in Cu02 planes, with infrared studies being inconsistent 
with the notion of static charge ordering [8| . 

We report an extensive study of the low-T (T~ 1 K 
and below) c-axis magnetotransport and low-frequency 
resistance (R) noise in LSCO with x = 0.03 [Fig. QJa)]. 
Such lightly doped samples are insulating at low T and, 
hence, most likely to exhibit charge glassiness 20]. The 
high-quality single crystal was grow n by the traveling- 
solvent floating-zone technique [2lJ; T sg ~ 7-8 K [9(. 
Two samples were cut along the main crystallographic 
axes and polished into 0.6 x 0.8 x 1.57 mm 3 (sample 1) 
and 0.6 x 0.9 x 1.6 mm 3 (sample 2) bars suitable for direct 
c-axis transport measurements. Both samples had the 
same behavior. R was measured with a standard four- 
probe ac method (~ 7 Hz) in the Ohmic regime, and noise 
with a five- probe ac bridge method [22l |. using a lock-in 
amplifier (~ 7 Hz; excitation current I exc = 1 nA in 
the Ohmic regime) to detect the difference voltage. This 
method minimizes the influence of T and I exc fluctuations 
on the R noise. Additional care was taken to ensure that 
the observed noise did not come from the contacts. In 
particular, the noise was measured twice on sample 2. For 
the first run, the contacts were made using the Dupont 
6838 Ag paste and, for the second run, by evaporating 
Au. In both cases, the room T contact resistances were 
< 1 fl, and the noise characteristics were the same. 

The c-axis R(T) exhibits insulating behavior, which 
obeys the variable-range hopping R = Ro exp(To/T) M for 
T < 1 K [Fig. W[h)]. The best fit to the data is obtained 
with fi = 1/3, consistent with early results on ceramic 
LSCO samples, where fj, = 1/2 for x = 0.02 and fi = 1/4 
for x = 0.05 [23j. Similar doping dependence of (i has 
been observed in various systems, such as doped semicon- 
ductors [3] and other disordered Mott insulators [25| . 
Surprisingly, here the precise form of R(T, B = 0) de- 
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FIG. 1: (color online) (a) Schematic phase diagram of hole- 
doped cuprates. (b) ZFC R(T). (c) Sample 2. The differences 
between field-cooled (FC) (squares: in 5 T, dots: in 9 T) and 
ZFC R vs. T vanish at Tb = T 5 x and TgT, respectively. 
Dashed lines guide the eye. The same results were obtained 
with 0.24, 0.014, and 0.007 K/min cooling rates. 




FIG. 2: (color online) Sample 1. R vs. time t upon the sub- 
sequent application and removal of (a) B — 1,2, 5, 6, 8 T, and 
(b) B = 6,1,2,8,1 T. (c) Sample 2; R vs. B. The arrows 
and numbers denote the direction and the order of B-sweeps. 
The sweep rate (0.001 T/min for B < 1 T, 0.005 T/min for 
B > 1 T) was low enough to avoid the sample heating. Inset: 
Subloops shifted vertically to for comparison. The hys- 
teretic and memory effects are similar for positive and nega- 
tive B, and for both B \\ c and B _L c. They become smaller 
with increasing T, and vanish for T > 1 K. 



pends on the cooling protocol. In particular, at low T, 
R(B = 0) obtained after cooling in field Bpc is higher 
than the zero-field cooled (ZFC) R(B = 0). This dif- 
ference decreases with increasing T, and vanishes at a 
temperature Tb that grows with Bpc [Fig-Ec)]- This 
history dependence, with nearly the same magnitude and 
Tb, was observed for both B || c and B _L c. 

Strong history dependent effects are seen also if B is 
applied after zero-field cooling. Figure [2^a) shows R as a 
function of time t upon the subsequent application and 
removal of several B values at a fixed, low T. Obviously, 
R increases with B. However, whenever B is turned 



off, R decreases, but it remains higher than the previ- 
ous R(B = 0). In fact, R(B = 0) keep increasing as 
long as subsequent B increase monotonically. Otherwise 
[Fig.[2[b)], R(B = 0) is determined by the highest previ- 
ous B: the sample acquires a memory of its magnetic his- 
tory. Similar memory effects in transport have been seen 
in other systems, such as manganites [26j in the regime of 
phase separation and heavily underdoped YE^CuaOg+x 
(YBCO) at low T [27]. In YBCO, they were attributed 
to the freezing of the directionally ordered charge stripes. 
In LSCO, the hysteretic and memory effects may result 
from the irreversible orienting of the FM moments [28| 
in the direction of the applied B. Therefore, it may be 
instructive to study the hysteresis in more detail. 

In order to obtain the data with a controlled history 
[e.g. Figs.|U[a)-(c)] the samples' memory was first erased 
by warming up above T~ 1 K (< T sg ), where hysteretic 
effects vanish [30j]. The subsequent cooling to the mea- 
surement T produced small R(t) relaxations only at the 
lowest T ~ 0.1 K, but any possible intrinsic nonequilib- 
rium dynamics could not be separated out from the ef- 
fects of cooling. The ZFC magnetoresistance (MR) R(B) 
[Fig. H[c)] first follows the paths 1-3 during the cycling 
of B between and 5 T. After the second sweep to 5 T, 
the first subloop (2-3) has closed, and continuing to raise 
B does not disturb the structure of the outer loop. The 
system thus exhibits return-point memory (RPM). The 
second closed subloop (5-6), obtained between the same 
B end points but with a different history, is clearly incon- 
gruent with the first one [Fig. WLc) inset] . This hysteresis 
in MR is strikingly reminiscent of the behavior of mag- 
netization in various magnetic materials [3lj ]. including 
spin glasses [32| . This suggests the existence of magnetic 
domains (clusters, or switching units) with holes confined 
within the domain walls, supported by the results below 
and consistent with earlier studies [1,0, H, HI, 13, 15. ljj]. 
The novel observations of the RPM and the incongruence 
of the subloops, which indicates that domains interact, 
impose strong constraints on theory. A model of an elec- 
tron nematic in Cu02 planes [33| finds the same type 
of hysteresis, but in the in-plane R anisotropy, so it is 
unclear if that can be related to our work. 

The noise, i.e. the time series of the relative changes 
in resistance AR(t)/(R), also provides evidence for the 
existence of some switching units. Figures |3Ja) and[3Jb) 
show clearly that, at low T, the system exhibits switching 
between well distinguishable states with the characteris- 
tic lifetimes ranging from a few minutes to several hours. 
Switching noise was also observed in YBCO, on much 
shorter t scales and at much higher T, in the pseudogap 
regime [34]. It is striking that here the noise amplitude 
is comparable to that in YBCO, even though our sam- 
ple volume is about 10 11 times larger. Also, here the 
switching noise is superposed on other random, slow fluc- 
tuations. The histograms of AR/(R) values [Fig. 12(c)] 
confirm that, at the lowest T (e.g. T ~ 0.1 K), the sys- 
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FIG. 3: (color online) Sample 1. (a) AR/(R) = (R-(R))/(R) 
({R} - the time-averaged R) vs. t for B — at several T. (b) 
An expanded section of (a). All traces are shifted for clarity, 
(c) PDF vs. AR/(R) for the 12-hour time interval, (d) The 
octave-averaged power spectra Sr(/) have been corrected for 
the white background noise. Solid lines are fits to Sr oc 1/ f a . 
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FIG. 4: (color online) Sample 1. (a) ZFC AR/(R) vs. t at 
T — 0.150 K for several B \\ c. All traces are shifted for clarity. 
Inset: A section of the noise at 5 T. (b) Sr(/ = 1 mHz) and 
(c) a vs. T (•) and B (□). The error bars in (b) are standard 
deviations of the data. Dashed lines in (c) guide the eye. 



tern appears to prefer two states, but other states are 
also present. Moreover, the precise shape of the proba- 
bility density function (PDF) depends randomly on the 
observation time (not shown), signaling that the system 
is nonergodic (glassy) on experimental time scales. All 
these effects become less prominent as T is raised, and 
the PDF becomes Gaussian already by 0.3 K. The noise 
power spectra Sr(/) (/ - frequency) obey Sr oc l/f a 



[Fig. OHd)]. ln order to compare the noise magnitudes 
at different T and £>, Sn(f = 1 mHz) is taken as the 
measure of noise, and it is determined from the fits 
in Fig. Eld). = 1 mHz) does not depend on T 

[Fig. 01b)], as i s already apparent from the variance of 
the raw data [Fig. [3£a)]. However, a increases from ~ 1 
to - 1.6 as T is reduced from 0.3 K to 0.1 K [Fig. gfc)], 
reflecting the slowing down of the dynamics and the in- 
creasing non-Gaussianity of the noise, similar to other 
systems out of equilibrium [ll| . 

The noise was measured also in B \\ c of up to 9 T. Sur- 
prisingly, B does not seem to have any effect on the noise. 
At low T, both switching events and other fluctuations 
on many different time scales are still present [Fig. IDJa)], 
the noise magnitude is independent of B [Fig. SJb)] , and 
Sr k 1//° with a ~ 1.6 that remains unchanged by B 
[Fig.|4fc)]. Furthermore, we have established that, unlike 
R itself (Figs. [T]and [2]), all the noise characteristics are 
independent of the magnetic history. 

The second spectrum S , 2(/2,/), which is the power 
spectrum of the fluctuations of Sr(/) with t, probes the 
correlations between fluctuators: it is white (independent 
of fa) for uncorrelated fluctuators (Gaussian noise), and 
Siifa, f) l//2~ /3 for interacting ones [13, [H, [36| . S2 
was calculated for a few octaves / = (/l, 2/^) [Fig.[5^a)]. 
The results [Fig.^b)] show clearly an increase of (1 — (3) 
from « at 0.3 K to large nonwhite values as T is re- 
duced, demonstrating that fluctuators become strongly 
correlated with decreasing T. On the other hand, B has 
no effect on large low-T values of (1 — (3), i.e. B does not 
seem to affect the nature of the correlations. 

As in studies of other glasses floL [ll| , we explore the 
scaling of S2{fa, f) with respect to fa and / in order to 
distinguish hierarchical pictures from generalized models 
of interacting, compact droplets or clusters. In the latter, 
the low-/ noise comes from a smaller number of large el- 
ements because they are slower, while the higher-/ noise 
comes from a larger number of smaller elements that are 
faster .10]. In the presence of short-range interactions, 
big elements are more likely to interact than small ones 
and, hence, non-Gaussian effects and S2 will be stronger 
for lower /. Indeed, we find that, unlike some well-known 



spin (e.g. CuMn) [lOj or Coulomb [llj glasses, S2(fa, f) 
is not scale invariant, but rather it decreases with / at 
constant fa/ f [Fig.[5jc)], consistent with droplet models 



37| . This result strongly supports the picture of spatial 
segregation of holes into locally ordered, interacting re- 
gions, as inferred from the MR hysteresis, resulting from 
competing interactions on different length scales. 

The noise study indicates that, at low T, the system 
wanders collectively between many metastable states. 
The apparent insensitivity of the noise statistics to B 
and to the magnetic history strongly suggests that the 
observed glassiness reflects the slow dynamics of charge, 
not spins. In conventional spin glasses, by contrast, all R 
noise characteristics are affected by B [38(. In addition, 
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FIG. 5: (color online) (a) The normalized second spectra 
<S l 2(/2), with the Gaussian background subtracted, for two 
T. Solid lines are fits to S 2 oc 1/ fl~ f) . (b) (1 - (3) vs. T (A; 
sample 1) and B (•: sample 1, ■: sample 2). The dashed line 
guides the eye. The results shown in (a) and (b) were obtained 
by averaging over the 0.5-1, 1-2, and 2-4 mHz octaves. The 
error bar in (b) shows the maximum standard deviation of 
the data, (c) 52 measured in octaves / = (/_l,2/_l). Dashed 
lines are linear fits to guide the eye. 



unlike spin glasses, here the onset of glassiness in both 
transport and noise occurs at T <C T sg , providing further 
evidence that the two phenomena are not directly related. 
In fact, a gradual enhancement of the correlated behavior 
with decreasing T suggests, in analogy with 2D Coulomb 
glasses [n}, that the charge glass transition T cg = 0. The 
noise reveals glassy dynamics only at T < 0.3 K, consis- 
tent with the higher T observations of static, as opposed 
to dynamic, charge ordering in underdoped cuprates in 
other experiments on comparable time scales [39( ■ 

We have demonstrated that the c-axis transport and R 
noise in x — 0.03 LSCO exhibit several clear signatures 
of glassiness, such as hysteresis, memory, and slow, cor- 
related dynamics, as T — > 0. The data strongly suggest 
that the doped holes form a cluster glass state as a result 
of Coulomb interactions albeit in the presence of a 
random, internal magnetic field that is produced by spin 
freezing. Further work at higher x is needed to determine 
how this dynamic charge order evolves into and coexists 
with the SC state. Studies of other cuprates, including 
those grown epitaxially, would be also of interest to see if 
the emergence of a charge glass is universal to cuprates, 
and perhaps even to all doped Mott insulators. 
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